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The scale factor a^s is the effective cross section used to characterize the measured 
rate of inclusive double dijet production in high energy hadron collisions. It is sensi- 
tive to the two-parton distributions in the hadronic projectile. In principle, the scale 
factor depends on the center of mass energy and on the minimal transverse energy 
-E'T.min of the jcts Contributing to the double dijet cross section. Here, we point out 
that proton-proton collisions at the LHC will provide for the first time experimental 
access to these scale dependences in a logarithmically wide, nominally perturbative 
kinematic range 10 GeV<ii'T,mm^lOO GeV. This constrains the dependence of two- 
parton distribution functions on parton momentum fractions and parton localization 
in impact parameter space. Novel information is to be expected about the trans- 
verse growth of hadronic distribution functions in the range of semi-hard Bjorken x 
(0.001<.x<0.1) and high resolution . We discuss to what extent one can disentan- 
gle different pictures of the x-evolution of two-parton distributions in the transverse 
plane by measuring double-hard scattering events at the LHC. 

I. INTRODUCTION 

In high-energy hadronic collisions, more than one pair of partons can interact with large 
momentum transfer. Such multiple hard interactions within the same hadronic collision become 
more numerous with increasing center of mass energy. They are a novel and generic feature 
of hadronic interactions at Tevatron and at the LHC. For instance, the inclusive cross section 
for double dijet production (see Fig. [Tl) is ctd (-ET.min = 20 GeV) ^ 10 /ib in proton-proton 
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FIG. 1: Schematic view of double dijet prodution in a proton-proton collision. 

collisions at the LHC, if each of the four jets carries more than a minimal transverse energy 
of -ET,min = 20 GeV. Even if this threshold is raised to -ET,mm = 100 GeV, the double hard 
scattering process is still in experimental reach with a^^^^ {E^^min = 100 GeV) ~ 50 pb. We 
will present calculations, supporting these estimates, in section [Tll 

It is not known how factorization theorems for large momentum transfer processes could be 
extended to multiple hard processes within the same hadronic collision, such as double dijet 
production. Typically, one assumes that such processes can be described as the incoherent 
superposition of single hard scattering processes 

mm 

. The double dijet cross section 
can then be expressed as the convolution of two independent hard partonic subprocesses with 
two-parton distribution functions Fd (see eq. (12.3p below). The ratio of the square of the dijet 
cross section as to the double dijet cross section for two indistinguishable hard processes 
defines the effective cross section aes 

^cfr = ^. (1.1) 

For two distinguishable hard processes A and B, it takes the form cXeff = ^^r^- If the two- 
parton distribution functions factorize into an uncorrelated product of standard single 
parton distribution functions, then cXefr gives access to the geometrical extension of the parton 
distributions in transverse space. 

There are several motivations for studying double hard cross sections at hadron colliders. 
First, hadronic collisions with more than one hard partonic scattering can contribute to multi- 
parton final states at high transverse momentum. Their improved understanding may help 
to control the QCD background to searches for novel physics in channels involving multiple 



high-£^r parton final states [8||, although simple kinematic cuts can be efficient in cleaning 
the signal from double dijet background sl- Second, multiple independent hard scatterings at 
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lower momentum transfers [Q^ ~ O ((1 — 5)^ GeV^)] play an important role in modeling the 



underlying event in hadronic collisions at collider energies [9|, llO, lUj. Studying the physics 



-n imp( 

0, m 



of such multiple hard scatterings at larger momentum transfer Q or as a function of Q may 



help to constrain the input to this modeling of the underlying event 



12| . Moreover, the double 



dijet cross section an provides qualitatively novel information about the transverse structure 
of the hadronic projectile because on depends on the relative transverse distribution of the two 
partons in the hadronic projectile 13|. If the two partons were distributed homogeneously over 
the entire hadronically active transverse area of the hadronic projectile, then the scale factor 
(Jeff would be comparable to the total inelastic cross section. The much smaller numerical value 
(Teff = 14.5 ± l-Tlll'l measured by the CDF Collaboration [3| disfavors such a homogeneous 
distribution. It is in support of a picture of the proton and anti-proton, in which partons with 
large momentum fraction are localized in a significantly smaller transverse region within the 
proton (see discussion of Fig. [2] below). 

Quantum chromodynamics offers a specific picture for the transverse growth of hadronic 
wave functions with increasing In 1/x or center of mass energy. In analogy to QED, where the 
transverse extension of the Weizsacker- Williams field of quasi-real photons around an electric 
charge grows with increasing energy, the hard (i.e. large-x) color charges in a QCD projectile 
can be viewed as sources of non-abelian Weizsacker- Williams fields, whose transverse size grows 
with increasing Inl/x or center of mass energy. The simplest, perturbative realization of this 
phenomenon in QCD is the BFKL evolution equation, which predicts with increasing \nl/x not 
only a growth of parton density locally in impact parameter, but also a growth of the hadronic 
projectile distribution function in impact parameter space, see e.g. Refs. 
Unlike QED, this perturbative picture is expected to be modified in QCD by saturation effects, 
which tame the growth of parton density locally in impact parameter, and by non-perturbative 
effects, which amputate the gluonic Weizsacker- Williams fields at a transverse distance set by 
confinement. One expects that this reduces eventually the growth of the average transverse 
extension of non-abelian Weizsacker- Williams fields from the perturbatively predicted power- 
law dependence oc x~'^ to a non-perturbative logarithmic increase. The scale and dynamics 
of the transition between perturbative and non-perturbative regime, as well as the physics in 
the non-perturbative regime remain under debate. All arguments indicate, however, that the 
transverse extension of the hadronic densities continues to grow at ultra-relativistic center of 
r.^. e„er,ie. albei. pos^W weake. than p.ed.cted pe.u.bat.vely U^U- 
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Here, we investigate to what extent this quahtative picture of the growth of transverse 
hadronic distributions could be tested by measuring double dijet production at Tevatron and at 
the LHC in the range of semi-hard momentum fractions, say 0.001<x<0.1, and hard momentum 
transfers. The starting point of our work is the observation that measurements of the inclusive 
double dijet cross section at the LHC can be performed over a wide range in -ET,min- In the 
following, we quantify this range and we investigate to what extent it provides access to a 
possible dependence of the scale factor on the cm-energy and transverse momentum cut-off. 
On the basis of these calculations, we assess the sensitivity of measurements of an to the 
transverse growth of the distribution of semi-hard partons in the hadronic projectile. 



II. THE FORMALISM 




We first discuss the formalism, on which the following calculations of double dijet production 
cross sections are based. To leading order (LO), the single scattering cross section to produce 
two massless partons of transverse energy larger than -ET,min reads 

cr siET, mm) = / dxi / dx2 / di^fi{xi,Q'^)fj{x2,Q'^)—^, (2.1) 

where 

' ' ^ 4 -^T min 

xi X2 > Xmin = and = p\. (2.2) 

The sum is over all possible quark and gluon 2— >2 scattering channels ij. For the purpose of 
the present study, expression f l2.ip is a sufficiently good approximation for the cross section 
of a dijet, produced in a single hard scattering event. We shall convolute LO parton-parton 
scattering cross sections with the CTEQ6L set of parton distribution functions [l^ for /j, 
which have been optimized for the use in LO calculations. The scale in the PDFs and in 
is set to the transverse momentum p\ = ut/ s for massless partons in the 2^2 process. We have 
checked that the results for cr5(-ET,min) obtained with this input from Eq. (12.1l) are consistent 
with Pythia 6.419 fl. 

We calculate the inclusive dijet production cross section as an incoherent superposition of 
two hard scattering processes within the same hadronic collision 

(Td (-ET.min) = 2 j '^^I'^^^dh j dxidx2dxsdx4diidi2 F}^{xi, X2]h — Si,h — S2) 

ijkl 
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xF^(x3,X4;si,S2) -J — 3^. (2.3) 

ati 0,12 

Here, the sum Ylijki '^^^^ parton species, which contribute to the two partonic processes 
i+j ^2 jets and k + l 2 jets. The symmetry factor | accounts for the fact that both partonic 
processes are indistinguishable in the sense that there is no operational prescription which 
establishes a one-to-one mapping between one dijet and one specific hard partonic cross section 
in f l2.3p . The Mandelstam variables ti and t2 are defined for two partonic 2^2 processes with 
incoming parton momentum fractions Xi, Xs and X2, X4, respectively. We choose the kinematic 
boundaries in the integrals over the incoming parton momentum fractions Xi, ... , X4 such that 
all outgoing partons carry more than a minimal transverse energy i?T,min. (In principle, more 
sophisticated kinematic boundaries could be implemented e.g. to require different i?T,min-values 
for both pairs of jets, but we shall not explore such possibilities in the following.) The integral 
in equation ( 12. 3p includes the two transverse positions Si, S2, at which the two hard processes 
take place, and the impact parameter b of the hadronic collision. The spatial information about 
the partons is specified in the two-parton distribution functions 

F},\xuX2-M,h2;QlQl), (2.4) 

which depend not only on the transverse momentum fractions Xi, X2 and the virtualities Qi, Q\ 
of both partons inside the hadron, but also on their transverse positions bi, b2. In the following, 
we shall often use a simplified notation, in which the virtualities are not written explicitly as 
arguments of Fq. For these virtualities, we will always choose the squared transverse momenta 
in the corresponding partonic 2^2 subprocess, as in (12. ip . (12. 2p . 



A. A factorized ansatz for two-parton distribution functions 

The discussion in this section is based on a class of models, which satisfy the factorized 
ansatz 

F'^{xr,X2; bi, bs; Ql Ql) = F\x^, bi, Ql) F\x2, ba, Ql) , (2.5) 

where 

F\x,h,Q')=nix,h)r{x,Q'). (2.6) 

A set of correlated two-parton distributions, which do not satisfy (12.5P will be discussed in 
section HVl In equation (]2.6I) . n{x,h) denotes the density of partons in the transverse plane. 
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It is normalized to unity, / dhn{x,h) = 1, so that f\x,Q'^) = J (ib b, Q^) are the stan- 
dard single parton distribution functions. For the class of models studied here, the transverse 
part of the density does not depend on the parton species i. The ansatz of Eq. (12. 6p con- 
tains information about the average transverse distance of the partons from the center of the 
proton in transverse space. If the partons are uncorrelated in impact parameter, the average 
distance between the two partons satisfies ((bi — b2)^) = (bf) + (b2). In this sense, the ansatz 
(12.51) also contains information about the average transverse distance between the two partons. 
In addition, it allows for a non-trivial x-dependence of the transverse size of the hadronic 
projectiles. 

Let us consider first the simple case of an x-independent density n{x, b) = ?T,(b). In this case, 
the geometrical information entering cross sections can be expressed in terms of the nucleon 
overlap function 

T^^ih) = j dsn{s)n(h-s) . (2.7) 

The normalization of n implies that / (ibTNN(b) = 1. For the factorized ansatz (12. 5p . the 
double dijet cross section then takes the form 

O'D (-f'T.min) = (^T,min)]' I J dhT^M, (2.8) 

and the scale factor ( 11.11) reads 

In the general x-dependent case, the nuclear overlap function (I2.7P will depend on the momen- 
tum fractions Xi,..., x^ of the partons in both hadrons. As a consequence, the scale factor 
becomes a function of the center of mass energy y/s and the jet energy threshold ii^T,mm- 



B. Interpretation of the scale factor in the model (I2.5p for two-parton distributions 

For the class of models (12. Sp . the scale factor contains information about the transverse 
parton density n{h) in the proton. A model-independent understanding of the b-dependence of 
this density distribution is missing so far. Here, we consider three-dimensional density profiles, 
from which transverse densities are obtained by projection, n(b) = f dzn{r). Rather than 
motivating a particular b-dependence of n(h), however, we prefer here to demonstrate that 
our conclusions about the scale factor will depend mainly on the rms of n(r) and will be 
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FIG. 2: The scale factor aes as a function of the 3-dimensional rms of the parton density distribution 
for a Gaussian, exponential and 6-profile n{r). The thinly dotted line denotes the central value of the 
Tevatron measurement Ces = 14.5 mb. The value (r^)^/^ = 0.875 fm of the proton charge radius is 
also indicated. 

rather insensitive to details of the functional shape. To establish this point, we compare three 
3-dimensional parton densities n(r). In particular, we consider a Gaussian density 
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(27r52)3/2 

a homogeneous density distribution which is sharply cut at radius R, 

3 



(2.10) 



and an exponential profile 
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87rA3 



exp 



|r| 

T 



f2.111 



f2.12l 



which has a more pronounced tail than the Gaussian distribution. To compare the sensitivity 
of (Jeff on the functional shape of the r-dependence of the density profiles f l2.10p . (12.111) and 
(I2.12p . we express the results in Figure [2] and Table [T] in terms of the root mean square (r^)^/^ 
of these densities. 

One sees in Figure[2]that the numerical value of (Jeff characterizes mainly (r^)^/^ g^j^^j ^j^g^^ j^g 
sensitivity to the detailed geometrical profile of n(r) is rather weak. As a consequence, the use 
of the Gaussian ansatz in the following studies can be regarded as a convenient choice which 
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Model for 3D density 


ms radius (r^) 




(r2)TEVATRON 


Gaussian oc exp(— ^) 


352 


3 ^ ' 


0.42 fm 


hard sphere oc Q(R — r ) 


5^ 


9.0 (r2) 


0.40 fm 


exponential oc exp(— ^) 


12A2 


7.3 (r2) 


0.45 fm 



TABLE I: Results from the spatial density analysis. Calculations of ms radius (r^) and aes are shown 
for several density models. The last column is calculated by equating the cJeff in the 4th column to the 
TEVATRON measurement of 14.5 mb. See also Fig. El 



will not bias our conclusions. We mention as an aside that a Gaussian profile may be motivated 



for instance by a study based on a light cone Hamiltonian 



2l| 



Figure [2] also demonstrates that the central value a^s = 14.5 mb of the CDF measurement can 
be related to a narrow range around (r^)^/^ ^ q 4 _ o.45fm for all three geometrical profiles. 
We note that if the partons relevant for double dijet production were distributed in impact 
parameter over a transverse region of the size of the proton's electric charge (r^)^/^ _ 0.875 fm, 
then the scale factor a^fi would take values between 50 mb and 70 mb, which are comparable 
to the total inelastic cross section (Xinei — 80 mb measured at Tevatron. A significant difference 
between the total hadronically active transverse size of the proton and the transverse extension 
of the region relevant for processes of high momentum transfer has been discussed repeatedly 



in the literature lla, \2 



27||. We note that the measurement of aes can not only provide an 



independent characterization of this difference. Moreover, it may also provide novel access to 
the dynamical origin of this difference via an analysis of the x-dependence of aes, to which we 
turn now. 
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C. Modeling the small-x evolution of the transverse size of hadronic wave functions 

The transverse size of hard partonic components of the proton is found to be smaller than 
total cross sections but will also grow with increasing ^/s. Therefore a larger part of the 
hadronically active regions in transverse space can be expected to contribute to hard production 
processes with increasing energy. 

We model this picture of the transverse growth of hadronic wave functions by specifying an 
x-dependence of the Gaussian density (12.111) . 

"(^'^) = 2^^"P 

We consider two parametrizations of the x-dependence of the Gaussian width 5{x). One model, 
originally proposed by Burkardt takes 

5{x) = wi ^{1 -x) ln(l/x), wi = 0.149 fm. (2.14) 

This results in a growth of the scale factor cTcs^x) oc 6'^{x) oc In 1/x oc In s, which is even weaker 
than the growth oc (Ins)^ which is realized in the Froissart bound. We fix the prefactor at 
wi = 0.149 fm to reproduce the Tevatron value for Ccs at -ET,min = 20 GeV. 

We also consider a second model which results in a power-law growth of aes oc s"^ with the 
center of mass energy, 

5{x) =W2 (l-x) x'"^/^, 1172 = 0.175 fm. (2.15) 

A power-law growth of the transverse hadronic wave function is obtained in the perturbative 
small- a; evolution, where the leading order BFKL-intercept is lj = ^ Nc 4: ln2. This is known to 
overestimate the growth of total cross sections in the experimentally accessible regime. However, 
it is conceivable that the growth of hard components in the transverse wave function is more 
rapid. Within the window of physically reasonable parameters, we have chosen u = 0.265 
to arrive at a model with power-like growth at small x. The prefactor in (I2.15P is fixed to 
W2 = 0.175 fm such that = 14.5 mb at i/i = 1.8 TeV and -ET.min = 20 GeV. 

In Fig. [3] we show the transverse growth of the width (b^) = 26'^{x) as a function of mo- 
mentum fraction x for the models (12. 14^ and ( 12.1511 . We note that a transverse growth similar 
to the model ( I2.14p has also been obtained in other calculations which model non-perturbative 
effects [18i l. 



26{x) 



(2.13) 
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FIG. 3: Two models for the increase of (b^) = 25^ (x) towards small x of the transverse proton wave 
function as a function of the parton momentum x. The lower (logarithmic) curve is a model (|2.14l) 



proposed by Burkardt 



22], the upper (power- like) curve is discussed in eq. (|2.15ll . 



A comment about our use of CDF measurements is needed here: in the CDF publication 
of 1993, a central value = 12.1 mb was quoted for -E'T.min = 18 GeV on the parton level. 
The later CDF measurement [7| quotes a central value a^g = 14.5 mb vv^ith much improved sta- 
tistical and systematic uncertainties. Hov^^ever, both CDF measurements considered a channel 
v\rith three jets and one photon in the final state, vv^hile v^e focus on four jet processes in the 
present discussion. The possibility that the scale factor differs for different channels due to 
their dependences on different parton distribution functions, cannot be excluded and has been 
explored in a recent model study [23|. To arrive at a particularly simple model, and since there 



is little phenomenological guidance, we do not explore the possibility of such differences here. 
We also remark that the more recent CDF analysis uses several lower values of -ET.min- For 
these reasons, we emphasize that our choice aes (v^ = 1.8 TeV, i?T,min = 20 GeV) = 14.5 mb 
does not reproduce the cuts and conditions of the CDF analysis. In particular, our choice of 
(Teflf = 14.5 mb can only be related to the CDF analysis under the assumption that the scale 
factor does not depend on the production channel. Moreover, our choice of -ET,min is dictated 
by the need of anchoring our discussion at a sufficiently large transverse energy ET,mim 

where 

the perturbative cross sections used in our calculations are sufficiently reliable. 
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III. NUMERICAL RESULTS FOR cJcfr FROM DOUBLE HARD CROSS SECTIONS 

In this section, we first characterize the range of i?T,min, which is experimentally accessible 
with sufficiently large event samples for the study of double dijet production at Tevatron and 
at the LHC. We then discuss the scale dependence of in the model (12.141) of Burkardt and 
in the BFKL model (12.151) . Finally, we turn to the question to what extent more complicated 
geometrical arrangements of two-parton distributions, or correlations not encoded for in the 
ansatz (12.50 can affect our conclusions. 

A. Rate of inclusive double 2-jet processes 

Figure [H shows the calculated cross section as for inclusive 2-jet production and an estimate 
of the double dijet cross section an as a function of the minimal transverse energy of the 
jets. To arrive at this estimate, the double dijet cross sections in Figure [4] is calculated 
from the single inclusive cross section as using equation ( 11.11) with a scale-independent value 
(Teff = 14.5 mb. If (Teff chaugcs with y/s, then the inclusive double 2-jet cross section an would 
differ from the value shown in Fig. [4] by a factor 14.5 mb/crcfr(-\A)- A table with numerical 
values for double dijet cross sections, entering Figures [4] and [5], is provided in the electronic 
supplement to this paper. 

At the Tevatron Run I (-y/i = 1.8 TeV), the inclusive double dijet cross section reaches 
^ 20 nb for -ET.min — 20 GeV. Upon increasing the jet energy threshold, this cross section drops 
rapidly to ^ 20 pb for -ET.min — 40 GeV. The kinematical reach at the LHC {^/s = 14 TeV) is 
much wider. With a cross section of ~ 16 nb, one gets to i?T,min — 50 GeV, and with a cross 
section of 10 pb, one explores the scale dependence of an up to -ET.min — 120 GeV. 

To put these cross sections into perspective, let us assume an integrated luminosity for Run 
II at Tevatron of 10fb~^ (more than 6fb~^ have been delivered to date). This would translate 
into 2 X 10^ double dijet events with -ET,min — 20 GeV and 2 x 10^ events with -ET.min — 40 
GeV. (For the purpose of these order of magnitude estimates, we have neglected the difference 
in center of mass energy between Tevatron Run I and Run II.) To relate the double hard cross 
section to a measurable quantity, it is necessary to disentangle the 4-jet events originating 
from two independent hard scattering processes from those that originate from a single hard 
scattering. There are various experimental handles for doing this. For instance, one can exploit 
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that in contributions to a^i, pairs of two jets must be balanced in Et, while contributions from 
other classes of 4-jet events are not and differ in the shape of their distribution. However, 
our study does not provide a basis for judging the (experiment-specific) size of event samples, 
needed for a measurement of a^- Here, we assume that ~ 10^ raw events are sufficient to this 
end, and we thus estimate that experiments at Tevatron Run II can measure aD for -ET.min ^ 40 
GeV. 

Under analogous assumptions, experiments at the LHC will access the physics of double 
dijet production for -ET.min ^ 120 GeV with the first 10fb~^ of data at y/s = 14TeV. 




50 100 150 50 100 150 



[GeV] [GeV] 

FIG. 4: Calculated inclusive single dijet cross section as (upper curves) and estimated double dijet 
cross section cjd = 173/(2 CTeff) (lower curves) as a function of the cut on jet transverse energy. Left 
hand side: as and ctd at ^/s = 1.8 TeV. Right hand side: as and ctd at -y/s = 14 TeV. The scale factor 
is assumed to be a^fi = 14.5 mb, independent of -ET.min and -y/s. 



B. The scale dependence of the scale factor cJeff 

The estimates given in section [Til Al above illustrate that experiments at the LHC can inves- 
tigate the scale dependence of the scale factor aes over a logarithmically wide range in £''p^min- 
We now illustrate the physical information which can be extracted from this scale dependence. 
To this end, we have calculated the inclusive double dijet cross section (I2.4p for two-parton 
distribution functions (12.51) with Gaussian transverse density profile n{x,h). Figure [5] shows 
results for the corresponding scale factor (Xcfr, calculated as a function of y/s and -ET.min for the 
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cases of a logarithmic and a power law x-dependence of the width S{x) of the parton density 
n{x, b) in the hadronic projectile. 




Et'"'" [GeV] Et^'" [GeV] 



FIG. 5: Theoretical calculation of the scale factor cjefr = (t| (ii^T.min, /2 cr_D (-ET.min, in the 
plane of center of mass energy and jet energy threshold £'T,min- Constant CTeff-values indicated on 
the right of the two plots lie on approximately straight lines. Inclusive single and double dijet cross 
sections are calculated from (12. 1|) and (|2.4I) respectively, for a Gaussian transverse density distribution 
of partons in the proton wave function. The x-dependence of the width of these density distributions 
is taken to follow a logarithmic increase (I2.14p [plot on left hand side] or a power-law increase (|2.15ll 
[plot on right and side]. Model parameters are chosen such that a^g = 14.5 mb for ^/s = 1.8 TeV and 
£'T,mm = 20 GeV. 

We expect that experimental data on a^s at the LHC will first become available as a function 
of -ET,min for fixcd ^/s. Howcvcr, the foreseen running schedule of LHC may also lead to 
information about the -y/s-dependence of double dijet cross sections. This is so, since before 
moving to = 14 TeV, LHC is scheduled to start operation this year with ^/s = 10 TeV, 
aiming for an integrated luminosity of 200 pb^^ which may be sufficient to explore double hard 
collisions over a range in i?T,min- Moreover, at a later stage in the LHC program, one may also 
expect a relatively short proton-proton run at ^/s = 5.5 TeV to collect comparison data for 
the LHC heavy ion program. In addition, data from Tevatron Run II may provide information 
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Ef"^ [GeV] 

FIG. 6: The scale factor cTeff for the two models of an x-dependent width of the Gaussian density 
distribution, shown in Fig. [3l Upper curves for ^/s = 14TeV, lower curves for ^/s = 1.8 TeV. 

about the -E'T.min dependence of an at y/s = 1.96 TeV. Measurements of the y^-dependence of 
(TD(-ET,min) will test the constancy of aes along lines of constant -ET,mm/ y/s- This constancy does 
not depend on details of the modeling of two-parton distribution functions, but results solely 
from the kinematic bound (12. 2p . This makes it an important consistency check for the picture of 
double dijet production advocated here. The main dynamical information of the measurement 
of (TefT is contained in its dependence on -Ex.min and y/s. The variation of this scale factor with 
kinematical variables reflects the growth of the transverse size of the projectile wave function 
with In 1/x. Fig. [5] illustrates this point for the two models of the transverse growth represented 
in Fig. [3l One sees that the scale factor increases significantly with increasing ^/s or decreasing 
-^T.mm- We emphasize that at y/s = 14 TeV, results of both models shown in Fig. [5] imply a 
variation of aes by roughly a factor 2 in the range between -ET.min = 10 GeV and -ET.min = 100 
GeV. This variation is much larger than the ~ 20% measurement uncertainty quoted by the 
CDF Collaboration for its measurement of aes- We take this as a strong indication that the 
x-evolution of the transverse size of hadronic wave functions in the range of semi-hard x is 
experimentally accessible via the measurement of double hard cross sections at the LHC. 

Once a non-trivial i?T,min-dependence of the scale factor a^s is established, one may ask the 
refined question of whether this allows for the discrimination between different models of small- 
X evolution. Since parton distributions f{x, Q"^) rise rapidly with increasing In 1/x, the inclusive 
single and double dijet cross sections (^^(ii^T.mm), <^ d{E-j: ,min) are dominated by x- values which 
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lie close to the lower bound of (I2.2p . As a consequence, the value of aes remains almost constant 
along lines of constant -ET.min/ in the (-ET.min, V^) plane. As an estimate, one may take 

aes^l5 6\x]ll), (3.1) 

where Xmin is the lower bound on the parton momentum fraction available for a dijet above 
-£'T,mm (see Eq. 12.21) ). Figure [6] shows in more detail that for the two models of x-evolution shown 
in Fig. m the steeper x-dependence in the transverse density profile n{x, b) is indeed reflected 
in a steeper ii^T.mm-dependence of cxefr. The model-dependent difference is more pronounced 
for the higher LHC center of mass energy and for lower values of i?T,min, where double hard 
scattering processes depend on parton distributions at smaller momentum fraction x. However, 
the differences are rather mild and may be difficult to disentangle experimentally. Moreover, the 
interpretation of relatively small variations in cXefr may require a more detailed understanding 
of the geometrical distributions entering (see section [TV] below). Thus, while a non-trivial 
-£'T,min-dependence of cXefr will allow one to disentangle models of small-x growth from the 
baseline of an x-independent transverse localization of partons in the proton, the discrimination 
between different models of small- a; growth may be more challenging. 



IV. CORRELATED TWO-PARTON DISTRIBUTIONS 



So far, we have discussed factorized two-parton distributions of the form (12. 5p . This ansatz is 
based on a picture, in which partons are centered around a single position b„ = in transverse 
space. The purpose of this section is to gain some understanding of the extent to which the 
results reached above depend on the geometrical assumptions underlying the ansatz (12.51) . and 
to what extent they refiect dynamical information. To this end, we shall study a simple model 
of two-parton distributions, which do not factorize into single-parton distributions. The model 
is based on picturing the transverse profile of the proton projectiles as being composed of three 
regions of size ~ 0.2 — 0.4 fm each, which have increased hadronic interaction probability. The 
centers b„- , i = 1,2,3, of these hot spots may be thought to be related to the positions of 



valence quar 
distribution 



is in the transverse plane. To be specific, we consider for these hot spots the 



25| 



(bti-^ , bjjj I 



exp 



(4.1) 
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Here due to the center of mass constraint the third coordinate is defined as b^g = — by^ — hy^ . 
We fix 6v = 0.25 fm, which corresponds to a separation of the centers of the three hot spots by 
an rms of 0.3 fm. Partons are located around the centers b„. of these hot spots at transverse 
positions b with Gaussian distributions 

= 2^ (-%:#)■ <^-2' 

We allow for an x-dependence of the Gaussian width 6s- A two-parton density distribution 
?T,/)(a;i, 3^2, bi, b2) for partons located at transverse positions bi, b2 can then be calculated in 
terms of the following integral over the centers of the hot spots b^^ , 

nz) (xi, 0:2; bi,b2) = ^ I (ib^irfb^al^lV'b^,^ 1^ ^(i(xi,bi,b^J(i(x2,b2,b^.). (4.3) 

Here, the sum over i and j averages over the different combinations of hot spots in one proton 
which can provide the two partons for the hard scattering process. In this class of models, the 
two-parton distribution functions do not factorize, but take the form 

F}^ (xi, X2; bi, bs; Ql Ql) = nD{xi,X2] bi, bs) f (xi, Ql) f (x2, Ql) . (4.4) 

The corresponding integration over transverse directions, entering (12.31) reads 

J dhdsi ds2nD (a;i,X2; b - Si, b - S2) ud ix3,X4; Si, S2) 

1 2 \ 



8n\6l 51 + 252 SI + 61 
where 

Si = dsixif + (5,(X2)' + 6s{xs)^ + 5,(X4)' . (4.6) 

In the limit 6^ ^ 0, one recovers the model of a density distribution with Gaussian profile, 
discussed in section UTTl More precisely, for 5^ — > 0, the centers of the three hot spots are all 
located at b„j = hy^ = b^g = 0, as can be seen from (14. ip . For the case of an x-independent 
density 5s(x) = 6s, Eq. (14. Sp becomes = g:^^ with (r^) = 3 5f. This is exactly the value 

of J dhT^j^{h) = l/o"efr for the Gaussian profile, obtained in Tab. [B 

To study the effect of non-factorizing distributions and compare to the results of section [1111 
we make use of the Gaussian width 5gjj = 6s{xY + 6^/3 of the corresponding one-particle distri- 
bution n(xi,6i) = / (ib2 n£)(xi, X2; bi, b2). We connect the correlated one-particle distribution 
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FIG. 7: Same as Fig. [5l but for a model of two-parton distribution functions, for which partons are 
localized in three hot spots in the parton wave function. 

of soft partons to the uncorrelated one by demanding that they have the same Gaussian width, 
6{x) = 6es- This yields 6s used for the numerical computation. 

We have calculated the inclusive double dijet cross section (I2.5p for the correlated two- 
parton distribution function (14.41) and compare again a power-law x-dependence of the form 
fl2.15p with a logarithmic x-dependence of the form (12.141) . The results of this calculation are 
shown in Figure [7l 

Similar to the calculation with factorized two-parton distribution functions (|2.5I) shown in 
Fig. [HI we observe from Fig.[7]that the scale factor acs increases with increasing ^/s or decreasing 
-E'T,min- The numerical differences between the results shown in Fig. [7] and Fig. [5] are relatively 
small. This may be understood by observing that the average distance between the centers of the 
ha*o„icany act.ve .e,o„. . ^^i^^Z^, ^ ^2 . 0.35 tn, f„. ^ 0.5 to. The ave.a.e 
distance from the center b^. , at which partons are localized is ^/(\^ = ^2 {5^/3 + 5^^{x)) > 
^^/2/35y. This is ~ 0.20 fm at large x and increases significantly for smaller x. As a consequence, 
the three hadronically active regions in the model MM overlap significantly for our choice of 
model parameters, and measurable properties of this distribution are likely to be similar to 
those of the single homogeneous density distribution, studied in section lllli 



18 



V. DISCUSSION 

In this paper, we have argued that the proton wave function is expected to grow in the 
transverse plane with increasing Inl/x for all values of x. We have then illustrated in model 
studies that this transverse growth should manifest itself in two-parton distribution functions 
and that it becomes experimentally accessible in the inclusive double dijet cross section ao- 
Of particular interest is the ratio (Xcfr = o's/'^'^d of the square of the inclusive single over the 
double dijet cross section. As a generic consequence of the picture advocated here, the scale 
factor (TefT is expected to be constant along lines of constant -ET.min/ V^, and it is expected to 
grow with a/s at fixed jet energy threshold -ET.min and to decrease with increasing -Ex.min at 
fixed y/s. We have shown that the wider kinematical reach of proton-proton collisions at the 
LHC will allow for the first time to test this scale dependence of the scale factor aes over a 
wide range in i?T,min- 

We have studied models for parton distributions in the transverse plane, which do not 
distinguish between gluons and valence and sea quarks of different flavor. In the absence of 
experimental information about such differences, we have adopted the baseline assumption that 
the geometrical distributions of all partons are the same, thus minimizing the number of model 
parameters. We note, however, that one may invoke QCD-inspired pictures for which different 
parton species are localized differently in the transverse wave function. In this case, the scale 
factor (Teff can depend on the production channel of the double hard scattering process, since 
different production channels (such as final states with two-jet events or with bb) depend on 
different parton densities. A model, which shows this feature, was studied for instance by Del 



Fabbro and Treleani in Ref. (23|. The study of specific double hard production channels has 
also been explored as a means to arrive at an improved (ideally: background free) experimental 
characterization of double parton collisions. In particular, the production of two equal sign W 



bosons at relatively low transverse momentum is dominated by double parton collisions [26|, 
though the cross section is very low. A much more abundant channel at the LHC, even if 
6-tagging efficiency is taken into account, is the production of bbbb, which may allow for an 
improved characterization of double hard cross sections [28 1. 



The models explored in the present paper result in an increase of with increasing ^/s. 



In contrast, the spatial correlations, implemented in the model of Ref. 



23|, imply that scale 



factors for all production channels decrease with increasing y/s. This illustrates that the models 
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studied in the present work do not exhaust all conceivable scale dependencies of cxefr. The 
models in section [1111 and [TV] are simple implementations of the picture that the transverse 
proton wave function grows with \nl/x. Moreover, the difference between our models and the 

n 

models of Ref. [23|| illustrates clearly, that far beyond producing only an experimental value for 
(Tefr, a measurement of inclusive double hard cross sections at the LHC can distinguish between 
qualitatively different pictures of the transverse proton wave function. As we have argued here, 
this may be done most efficiently by studying the -Ex.min-dependence of the scale factor aes at 
fixed y/s. 
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